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The heat kernel method is extended to the case of finite temperature. Special emphasis is given 
to the study of gauge theories. Due to the compactness of space in the Euclidean time direction 
(inverse temperature) the field strength cannot completely characterize the gauge fields. This is just 
a manifestation of the Aharonov-Bohm effect. The field strength has to be supplemented by the 
Polyakov loop. Only if the latter is taken into account one obtains gauge covariant results for the 
generalized Seeley-DeWitt coefficients of the heat kernel expansion. 
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O ! I. INTRODUCTION 

CM . 

. Finite temperature quantum field theory has gained widespread applications in different areas of physics like e.g. for 
r ^ ' the description of the early universe or the study of heavy- ion collisions. This work is devoted to a finite temperature 
generalization of a non-perturbative technique which is very well established for the vacuum case. The heat kernel 
OO ■ expansion [HQQ (for a review see e.g. 0) is frequently used in the context of the path integral method to integrate 
\ out degrees of freedom in a non-perturbative way. The outcome is a derivative expansion of the fields which correspond 
to the remaining degrees of freedom which were not integrated out. In other words the heat kernel expansion yields 
an effective field theory (cf. e.g. 0, El El)- A typical application which we also address here is to integrate out 
, fermions from the action of a gauge theory which is bilinear in the fermion fields. The obtained fermionic functional 
determinant is processed further by the heat kernel expansion. In this way the ultraviolet divergent terms which 
t-H , require renormalization can be properly isolated e.g. using the function method Q. It may appear that part of the 
symmetries of the original action are spoiled by the need for renormalization. Such an appearance of anomalies can 
also be well accounted for with the heat kernel method. In the following we shall not consider anomalies explicitly. 
We keep the number of dimensions as well as the gauge group completely general. We restrict ourselves, however, to 
a flat space-time structure. 

Naively generalizing the heat kernel expansion to finite temperature we will encounter problems concerning the 
gauge symmetries as we will see below. This, however, is not a manifestation of an anomaly, as it appears for any 
Oh, gauge symmetry and only at finite temperature. We will solve that problem subsequently by presenting an improved 
heat kernel expansion which preserves the gauge symmetries (apart from potential anomalies of course). Throughout 
this work we will use the Matsubara formalism 0, 13 to represent fields at finite temperature. It amounts to a 
Euclidean metric with the boson (fermion) fields obeying (anti-)periodic boundary conditions in the Euclidean time 
direction with periodicity (3 = 1/T where T is the temperature of the heat bath the fields are coupled to. 

The task to extend the heat kernel expansion to finite temperature has already been addressed in 9]. It has been 
claimed there that the ratio of the finite temperature heat kernel and the vacuum one is just a temperature dependent 
constant independent of the considered fields. In contradiction to that statement the authors of ^} presented a 
counter-example where the mentioned ratio does depend on the gauge fields. It is the purpose of the present work to 
clarify this contradiction and to present a valid generalization of the heat kernel expansion to the finite temperature 
case. In vacuum the heat kernel expansion is an expansion in powers of the field amplitudes and their derivatives. For 
gauge fields this amounts to an expansion in powers of the field strengths and their gauge covariant derivatives. All the 
other information present in the vector potentials is due to the gauge freedom and therefore of no physical significance. 
At finite temperature, however, the field strengths and their derivatives are no longer sufficient to characterize the 
gauge fields. The crucial point is that due to the periodicity requirement mentioned above the Euclidean time direction 
is compactified. This leads to the appearance of an Aharonov-Bohm effect The Polyakov loop 0] supplements 
the information contained in the field strengths. It is defined by 
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L(x) := Pexp 



Xa+P 

dtA (t,x) 



(1) 



where Aq is the zeroth component of the gauge potential and P denotes path ordering. The latter is required for 
non-abelian gauge theories to obtain a gauge covariant result. Obviously the Polyakov loop constitutes a non-local 
object in terms of the gauge fields. It is easy to check that a Taylor expansion of L in powers of the gauge potential and 
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its derivatives e.g. around the point x would yield objects which are ill-defined with respect to gauge transformations. 
Therefore, the Polyakov loop should be treated on equal footing with the field strength without any further expansion 
in terms of the amplitudes of the gauge potential. On the other hand, a heat kernel expansion is just an expansion in 
terms of the field amplitudes. Therefore to prevent an erroneous expansion of the Polyakov loop, the latter has to be 
isolated right from the start to obtain a heat kernel expansion which is appropriate for the finite temperature case. 

The paper is organized in the following way: In the next section we take the usual method to perform a heat 
kernel expansion and naively generalize it to the finite temperature case. To some extent we follow here the work of 
0. At a certain stage, however, we have to point out differences between our findings and the ones of Finally 
we have to face a problem concerning the gauge covariance of the results obtained with this naive generalization 
of the heat kernel expansion technique. In Section lllll we present an improved method to perform the heat kernel 
expansion at finite temperature preserving gauge covariance in every single step. Gauge covariant Seeley-DeWitt 
coefficients are calculated in Section llVl It is also shown in this section that the results reduce to the already known 
results for the vacuum case. In addition, the general results for the Seeley-DeWitt coefficients are applied to the 
case of two-dimensional QED with constant fields. For this special case the heat kernel and thus also all its Seeley- 
DeWitt coefficients can be calculated exactly ^3]- This provides the possibility to cross-check our results with the 
ones presented in [10J. Finally Section Ivl summarizes our findings. Three Appendices are added containing technical 
aspects of the presented formalism. 



II. NAIVE REPRESENTATION FOR THE HEAT KERNEL 



Consider the differential operator in Euclidean space 

A, = -D-D; + X(x) (2) 

with the gauge covariant derivative 

D^ = d; + A^x). (3) 

Here the functions and X are no operators (but might be matrix- valued due to the group structure of the gauge 
group under consideration). Under gauge transformations the introduced quantities transform as 

X(x)^V(x)X(x)V- 1 (x), — ► V(x) D^ l V~ 1 (x) . (4) 

The heat kernel is given by 

H(x, x'; t) := (x\e~ rA \x') = e' T ^ S(x - x') . (5) 

We shall assume that the differential operator J5J originates from a fermionic determinant which we want to evaluate 
by the heat kernel technique 0. A generalization of the presented formalism to the bosonic case is in principle 
straightforward. The representation for the heat kernel on the right hand side of (jSJ) is completely appropriate for 
the vacuum case. In this section we will generalize this heat kernel to the finite temperature case. We will point out 
the problem which appears when this heat kernel is used at finite temperature. In the next section we will present 
an improved representation for the heat kernel which is better suited for the treatment at finite temperature. Hence 
some of the formulae which we present in this section should not be regarded as valuable results but as illustrations 
of the problems which appear at finite temperature. 

For finite T the 5-function which appears in (|SJ is given by 

+°° r riD-ii, r rl D h 

S(x - x') = T j (f^T ^ S } ^ {X °~ Xo) =■ f (f^p ^' (X ~ X ] (6) 

where D is the number of space-time dimensions. Being interested in the calculation of a fermionic determinant via 
the heat kernel method the discrete energies appearing in JSJ are defined as Q, @ 

uj n = 2irT (n + 1 ) . (7) 



(8) 

/=i 



2, 

It is straightforward to show that the diagonal entries of H can be written as [l3T | 

H(x, x-r)=^ ±± a'* ^+ M ^ £ ± r m (2zfc ■ D x - A x + M 2 ) m f(x) 



m=0 



m 
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We have introduced an arbitrary mass term M to keep the formalism as general as possible. All derivatives in (jSJ 
act on an arbitrary function / which has to be put to 1 at the very end. We want to emphasize already here that 
the precise meaning of "covariant derivative acting on /" is important. As / basically substitutes a fermionic field 
the meaning of D^f(x) really is d^f(x) + A fJi (x) f(x) and not [D*, f(x)]. Therefore eventually putting / to 1, a term 
like Apfa) might remain! This is in contrast to the claim in |9j that D^f -terms should be dropped in the end (cf. the 
statement between eqs. (3.17) and (3.18) in Q). This will have important consequences as we will see below. 
Next we introduce generalized Seeley-DeWitt coefficients as follows: 
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D/2 



H & = ( 4— J e L 9n(^ V^T) r" . (9) 

^ ' n=0 

Note that at finite temperature there is a remaining r dependence of the coefficients g n but only via the specific 
combination \frT 9]. Hence the g„'s become independent of r for the vacuum case T = 0. To calculate e.g. the first 
two coefficients, <?o and g±, we have to evaluate the sum over m in (|8J) up to m = 2. We note that all terms with odd 
numbers of drop out by the k summation/integration. We first calculate 

c™ := (2ik -D x - A x + M 2 ) m f(x)\ f , ^ . (10) 

J ~ L drop odd powers of fc 

We get 

c = 1, (11a) 
ci = DpAn-X, (lib) 
c 2 = —Ak^kuD^Au + DuD^D^A^ - D^D^X - XD^A^ + X 2 , (11c) 

where we have introduced 

X := X — M 2 . (12) 

Note that the momenta scale as 1/y/r due to the exp(— rk 2 ) term in ©. Therefore it is easy to see that cq contributes 
to go, c\ and the first term of ci contribute to g±, and the remaining terms of ci contribute to 32- Hence 

.go = (4vrr)^/ 2 ^ ^ c = h{^T) , (13a) 

9l = M D/2 ^^je- rf U-^TAk^KD^ =(D fl A li -X)I (V^T)-D A I 2t (V?T)-D i A i I 2s (V7T) 

= -XIo(V^T) + D A (io(V^T) - ht{^T)) , (13b) 

where the relevant sum-integrals are collected in Appendix lAl 

The important point now is that the term DqAq does not transform homogeneously under gauge transformations in 
contrast to X. 1 Therefore this cannot be a proper result for the generalized Seeley-DeWitt coefficient. This problem 
does not show up at zero temperature since then both integrals I and I 2 t become identical (cf. Appendix Hence 
the dubious term drops out in this case. We note that this Do^o-term was overlooked in It was still present in 
eq. (3.15) of but was disregarded afterwards by arguing that D f vanishes in the end. As already pointed out, 
this would be true for [D ,f] but not for D f. Ignoring this subtle difference the authors of 9] did not encounter 
any problem concerning gauge invariance. However, as shown above there is a problem with gauge invariance and 
actually there should be one on physical grounds: By construction the heat kernel expansion JHJ is an expansion in 
powers of field amplitudes and their derivatives. Any non-local object is expanded in terms of local ones. In vacuum 
this is a procedure which preserves gauge invariance — apart from well-defined anomalies which should be there on 
physical grounds. The reason why this procedure works for the vacuum case can be traced back to the fact that the 
field strengths (and their derivatives) are sufficient to characterize the physical content of the gauge fields. The field 
strengths, however, are local objects. The situation is completely different at finite temperature as already mentioned 
in the introduction. In this case the Polyakov loop carries additional information. This object is manifestly non-local. 
An expansion of the Polyakov loop in powers of the field amplitudes would yield misleading gauge dependent results. 



Note that X transforms in the same way as X. 
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Hence the Polyakov loop has to be isolated as an additional degree of freedom before processing the heat kernel 
expansion. 

Having established the physical reason for the erroneous result (|13b|l we next have to analyze what went wrong on 
the technical level when proceeding from © to JSJ. Concerning the behavior with respect to gauge transformations 
we observe that formally the heat kernel (JSJ transforms as a non-local object, 

H(x,x';t) -> V(x) e~ TA * V^{x) S(x - x') = V(x) e~ T ^ 5{x - x') V- X {x') = V(x) H(x, x'; r) V' 1 ^') , (14) 

on account of the transformation properties I0J. It is illustrative to decompose this transformation property in the 
following way: 

e -rA. ->v(x)e- T *"V-\x) and S(x - x') -> V{x) 5{x - x') V~\x') . (15) 

Correspondingly the left hand side of (jHJ) transforms gauge covariantly 

H(x,x;t) -> V(x)H(x,x;t)V~ 1 (x). (16) 

In contrast, the single Fourier modes on the right hand side of iJHJ do not show that behavior with respect to gauge 
transformations. In general, only after summing up all Fourier modes and performing the sum over m one obtains 
the heat kernel with its well-defined transformation properties. For the vacuum case it turns out that already the 
summation over the Fourier modes is sufficient; each single Seeley-DeWitt coefficient obtained in that way transforms 
gauge covariantly. As we have demonstrated above this is no longer true at finite temperature. Therefore it is more 
appropriate to use an approach which preserves in each step as much of the good transformation properties as possible. 
To be specific, we are looking for a Fourier type representation of the 5-function which preserves the transformation 
property (|15J) for each single Fourier mode. In this way we will never encounter the problem that only a sum of objects 
transforms properly while the single terms do not. We stress once more that such an approach is not mandatory for 
the vacuum case. After the summation of the Fourier modes is performed the results transform properly with respect 
to gauge transformations. 

III. GAUGE COVARIANT HEAT KERNEL 

To obtain gauge covariant results for the Seeley-DeWitt coefficients it might be not mandatory but it is obviously 
more save to start with a Fourier type representation of the <5-function which transforms gau ge c ovariantly in every 
single Fourier mode. An important ingredient for such an approach is the gauge link (see e.g. |14lll5lll6l|) 

U[A](x;x') := Pexp 

where P denotes path ordering. For our purposes it is useful to choose the following path parametrized by z: the 
path connects with straight lines the sequence of points x' = (x' , x') — > (to, x') — > (to, x) — > x — (xo, x). Here To is an 
arbitrary time. It will turn out that none of the final results will depend on To. From the technical point of view this 
aspect can be seen as a consistency check of the whole approach. For the chosen path we find 

U[A](x,x) = l. (18) 

Another important property of gauge links is their additivity |l4l Il5l Il6j : 

U[A](x;x')U[A](x';x") = U[A](x;x") (19) 

valid for arbitrary x' on the contour connecting x and x" . Finally under gauge transformations U transforms as 

Guam 

U[A](x;x') -» V(x)U[A]{x-x')V-\x'). (20) 

Therefore it is tempting to write 

6(x - x 1 ) = S(x - x')U[A](x,x') = ^ ' } U[A](x, x') . (21) 



- dz^A^z) 



(17) 



5 



Indeed, now the modified single Fourier modes transform gauge covariantly, i.e. in the same way as the (^-function in 

e *<"-*')U[A](x,a/) V(x) e ik < x - x "> U[A](x;x')V- x {x') . (22) 

However this cannot be the full story: At finite temperature the 5-function which corresponds to fermions must be 
antiperiodic in Euclidean time with periodicity (3 = X/T: 

5(x +/3-x' ) = -5(x -x' ). (23) 

This is also true for the single Fourier modes 

gifco (x +/3-x' ) _ _gifco (x -x' ) ^24) 

on account of JJJ. The modified Fourier modes displayed on the left hand side of (|22J) do not share that important 
feature as the gauge link U is not periodic in time: 

U[A](x + P,x;xq,x') = U[A](xo + (3,x;x 0: x) U[A](xo,x]x' ,x) ^ U[A](xo,x;x' ,x) (25) 
s , ' 

=L(x) 

where we have used (|19|l . Note that the term which spoils the periodicity is nothing but the Polyakov loop defined in 
©■ 

We now construct a modified gauge link which both has the transformation property l|2(J|) and is periodic in time. 
The key point is to isolate the (generalized) zero mode of the time component of the gauge potential which is closely 
related to the Polyakov loop. We first define the generalized zero mode 

oq(x) := -TlogL(x) = -Tlog U[A](x + /?, x; x , x) . (26) 

The properties of ao are collected in Appendix [BJ In general already the log- function of complex numbers is multiple 
valued. This is of course also true for the log-function of matrices. In the following we use the log-function defined 
via its Taylor expansion (cf. Appendix [Bj 

log(l + 0) :=^2(-l) n+1 -O n . (27) 

71=1 71 

Hence we restrict ourselves to Polyakov loops with values close to the unity matrix. We postpone the further discussion 
of that issue to Section We introduce the modified gauge link 

U[A][x;x') :^ e Xoa "^ U[A](x;x') e- x 'o a ^ x "> . (28) 

Using the properties of ao collected in Appendix iBl it is straightforward to prove the following properties of U : 

U[A](x;x) = 1, (29a) 
U[A](x a + /3,x;x' ,x') = U[A}(x ,x;x' ,x') , (29b) 
U[A]{x;x') -> V(x)U[A]{x;x')V- 1 (x') . (29c) 

With these developments at hand the heat kernel can be represented as 

H(x, x'; t) = e- TA * (S(x - x') U[A]{x; x')\ . (30) 

For the diagonal part one gets 

H(x, x;t)=^ e- T ^ +M ^ ^ — r m (2ik ■ D x - A x + M 2 ) m U[A] (x; x') ( ?,\ ) 



ml 

m— 
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IV. GAUGE COVARIANT SEELEY-DEWITT COEFFICIENTS 



From the representation l|31|l it is straightforward to calculate the generalized Scclcy-DcWitt coefficients defined in 
©. We present in the following the calculation of the first three generalized Seeley-DeWitt coefficients. For this task 
we introduce 



(2ik ■ D x — A x + M 2 ) m U[A](x; x 1 ) 



drop odd powers of k 



(32) 



Note that by construction every c m transforms gauge covariantly in contrast to c m as defined in (|f (Jfl . For the lowest 
order terms we get 



co = 


1, 










Cl = 


(D^D^-X^A^x') 


1 

X —X 








C2 = 


{-Ak^k v D^D v + D^D^DyDy - XD^D^ - 


- D^D^X + X 2 )U[A]{x-x') 


x' — X 


C3 = 


Ak^K (-D^D„D a D a - D^D a D a D v - A 


D a D^D, + D^D„X + D^XD, + XD^D U ) U[A] (x; x 




+terms which contribute to <?3 , 








c 4 = 


16k fi k u k a k /3 D fl D u D a D U[A}(x] x') 


4 


terms which contribute to 53, (74 . 



(33a) 
(33b) 

(33c) 

—X 

(33d) 
(33e) 



It is lengthy but straightforward to determine the corresponding generalized Seeley-DeWitt coefficients. In the course 
of the calculation there appears the following type of n-fold derivatives: 



Dl U[A)(x;x') 



a%(x). 



(34) 



The proof for this relation is given in Appendix [CJ The formulae collected in Appendix ^ serve to perform the ki 
integrations and to rewrite the finite temperature sums. The lowest order generalized Seeley-DeWitt coefficients turn 
out to be 



.go - hWrT), 

9l = -X Io(V^T) + a 2 (lo{V^T) - I 2t (V^T)) , 



(35a) 
(35b) 



1 



f 



I r 



52 = -X 2 I a {^?T) - - [A, [A,*] J Jo(v^T) ~ g [D , [A,*] J (3J (y/rT) - 2/ 2t (V?T)) 



D , X] ao + Xaf) (lo(V^) ~ htiV^T)) + D,} 2 I Q (V^T) + ~[A, A] 2 (2J (V?T) - ht(V^T)) 



i [A, [D , [A, -Do]]] (ioCv^T) ~ litWrT)) + 5 [A, [A, A]] «o (/o(v^T) - MV^T)) 



(35c) 



This constitutes the main result of the present work. Of course, by construction every coefficient transforms gauge 
covariantly. Obviously the method can be extended to calculate higher order coefficients. The expressions for I are 
collected in Appendix [S] Note that the coefficients g n do not depend on the number of space-time dimensions D. 
The dependence of g n on the temperature is non-analytic as can be seen by inspecting (|A6|) and l|A7|) . This leads to 
a very smooth transition of the vacuum to finite temperature results. In the following we shall apply these general 
results to two special cases. The purpose is to check the consistency of our results. 

In vacuum all / -functions become unity (see Appendix lA^I . One obtains the standard results for this limiting case 
1: 



So™ 


= i, 






= -x, 




3 2 vac 


- 


I - 

6 . 
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Di,[Di,X} 



Do,[D ,X]\ + A[A,Af + ^[A,A] : 



1 



D„[D„X\\ +-[D„,D U ] 



(36a) 
(36b) 



(36c) 
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As a further check of our results (|35|l we compare them to the special case discussed in two-dimensional 
massless QED with a constant electric field strength E and in addition a constant zero mode. The (trace of the) heat 
kernel for such a system is calculated exactly in [lCj : 



txH(x, x; t) 



E 



1 



2-7T tanh Et 



l + 2^T(-l) n cos 



Ex\ + 2irTa 



T 



exp 



n 2 E 



AT 2 tanh Et 



(37) 



Below we will expand the exact result to obtain the corresponding generalized Seeley-DeWitt coefficients. The latter 
can be compared to our results which we will apply now to the case at hand. 
The vector potentials are given by 



A a = -i (Exi + 2nTa) , A x = . 
For the zero mode defined in (|26|l we find 2 

ao = Aq — —i (Exi + 2nTa) . 
The task is to calculate (approximately) the heat kernel for the operator 



fl 2 = D^D, 



X 



with 



X = --[^, 7u ][D fl ,D„]=-Ea 3 



where we have used the convention [Toj 



7o = — 01 , 7i = a '2 

with the Pauli matrices Cj. 

For this case the generalized Seeley-DeWitt coefficients are (of course we use M = 0): 

Ea 3 Io(V^T) - [Exi + 2iTTa) 2 (l {V^T) - htiV^T)) , 

-Ea 3 {E Xl + 2nTa) 2 {h{yf?T) - ht(^T)) + ^E 2 (I (t/tT) + ht(^T)) 
1 



9o 
91 

.92 



(38) 
(39) 
(40) 

(41) 
(42) 



(43a) 
(43b) 



(43c) 



The expansion of 13 7|) has to be performed in the following way: Expand in powers of r, but with T 2 t fixed. There 
is, however, one exception to this rule, namely the 2irTa term. The temperature appearing there docs not originate 
from integrating out the thermal fermions but was put in by hand in |l0| to make a dimensionless. To account for 
that fact we use (|39|l and rewrite (|37|l as 



tvH(x, x; t) 



1 



Et 



2-kt tanh Et 
Now we apply our expansion rule to l|44|) : 

txH (x, x; t) 



l + 2^(-l) n cos 



r- iao 1 



exp 



Et 



AtT 2 tanh^r 



(44) 



1 

2ttt 



1 + -t 2 E 2 



x < 
1 

2ttt 
1 

2ttt 



l + 2^(-l)"(l + r 

n=l 
1 

3 

1 



1 + -t 2 E 2 



2tT 2 24t 2 T 4 
1 



2 T 4 J V 3 4tT 2 J t V 4rT 2 J 



6 
1 



ra 2 - -t z E z (/„ - ht) + T:T Z al (J - 27 2 



1 



1 + -t 2 E 2 ) h + ( ra 2 - -t 2 E 2 ) (J - I 2t ) + -T 2 a% (J - 2I 24 + J 4i ) 



(45) 



2 The already mentioned ambiguity concerning the logarithm will be discussed below. 
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where we have used the relations (|A6|) and i|A7(l to rewrite the finite temperature sums. Finally recalling JSJ we can 
read off the (trace of the) generalized Seeley-DeWitt coefficients: 

tr 5o = 2/ , (46a) 

tvgx = 2a 2 (I -ht), (46b) 

tr 32 = l i E 2 {h+ht)+4{h-2I 2t + h t ). (46c) 
These results are in complete agreement with l|43|) ■ 



V. SUMMARY 



We have presented a generalization of the heat kernel expansion technique to finite temperature. This task has 
already been envisaged in Technically we followed the approach of to some extent. There are deviations, 
however, concerning the formal as well as the physics point of view - which are of course intertwined. Physically we 
expect to find (at least) two aspects which make the finite temperature case distinct from the vacuum case. First, 
we expect the appearance of terms which "break" Lorentz invariance like e.g. a different coefficient in front of the 
electric field strength (squared) as compared to the magnetic field strength. Such differences have not been found in 
whereas our calculation shows that they are indeed present (cf. e.g. the fifth and sixth term in (|35c(0 . Second, as 
already pointed out in the introduction at finite temperature the field strengths are not sufficient to characterize the 
gauge fields. Due to the compactness of the Euclidean time direction the Polyakov loop carries additional information 
— a manifestation of the Aharonov-Bohm effect. Hence we expect the appearance of a corresponding object in the 
heat kernel expansion. Indeed in our approach the generalized zero mode which is the logarithm of the Polyakov loop 
appears. Again this object does not show up in the approach of 9]. 

Concerning the technical aspects of our approach as compared to [t| we have taken care of the correct meaning 
of the gauge covariant derivative acting eventually on the arbitrary function /. The latter appears in Nepomechie's 
constructive method to obtain the heat kernel expansion; cf. JHJ. Only the correct treatment of the covariant derivative 
revealed the problem of the naively calculated Seeley-DeWitt coefficients concerning gauge covariance. Subsequently 
we developed an expansion method which preserves gauge covariance in every single step. While this is not necessary 
for the vacuum case — albeit it does not hurt — it turned out to be very useful for the finite temperature case. 
The requirement to simultaneously preserve gauge covariance and (anti-)periodicity led to the introduction of the 
generalized zero mode. In this way the Polyakov loop entered our formalism. 

Formally the results of could be reproduced if we replaced lit and I^t by I in . Obviously a large number 
of terms would then drop out. These terms have not been considered in 0- By inspecting l|A6(l it is clear that such a 
replacement is only correct at zero temperature. Consequently the terms which are neglected in |9j are the ones which 
are not Lorentz invariant. However, the appearance of such terms is expected as the heat bath "breaks" Lorentz 
invariance. On the other hand, rotational invariance is still preserved. Obviously this requirement is fulfilled by our 
result JSSJl. 

Further support for the soundness of our results comes from where the case of two-dimensional QED at finite 
temperature with constant fields is addressed. The possible appearance of a zero mode which cannot be gauged away 
at finite temperature constitutes the discrepancy between [j| and |Tol | as noted in the latter reference. As shown above 
our results are in complete agreement with the special case of two-dimensional QED with constant fields discussed in 
|10| . One might regard our results as a generalization of the results of |10| to arbitrary dimensions and non-constant 
fields. On the other hand, in contrast to |lOj| we cannot offer an exact solution here but only an expansion in terms 
of field amplitudes and their derivatives. 

Concerning gauge fields at finite temperature the crucial point is that this expansion is not only an expansion in 
powers of the field strength but in addition in powers of the generalized zero mode. As already emphasized in [Toj 
even a constant zero mode cannot be gauged away at finite temperature. It is important to stress one aspect which 
is different in our approach as compared to the exact result presented in 0] . Inspecting (|57J) reveals a periodicity of 
the heat kernel with respect to 

a -> a + 1 . (47) 

Our approximate result l|43(l seems not to reproduce that behavior of the exact solution. Of course, an expansion in 
powers of the zero mode amplitude might not be appropriate any more for large amplitudes. On the other hand, due 
to the periodicity property l|47|l the zero mode can be made arbitrarily large. Actually this problem is just another 
manifestation of the fact that the logarithm introduced in (|26[) is in principle multiple valued. Indeed the Polyakov 
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loop which corresponds to the problem studied in |Toj is given by 

L{x) = e *(/3^i+ 2 ™) . (48) 

Obviously at this stage the periodicity property is still present. The problem appears with the calculation of the zero 
mode according to l|26() which apparently yields l|39|) . This, however, is not fully true: If the logarithm is defined via 
its Taylor expansion (|27fl the outcome of l|26|l is 

a (x) = -i [Exi + 2nT(a - m)] (49) 

where m is an integer chosen such that <f> := \(3Ex± + 2ir(a — m)\ is minimal. And even this is only true if <j> is 
small enough, i.e. |e^ — 1| < 1. Otherwise the Taylor expansion (|27J) for the logarithm would not converge. With 
the correct result (14911 for the zero mode we see that the periodicity property (|47|l is still present. On the other 
hand, by considering the convergence properties of the Taylor expansion for the logarithm it becomes obvious that 
for an expansion of the heat kernel small amplitudes are mandatory right from the start. Strictly speaking so far 
our formalism can only be applied to Polyakov loops with values close to unity. Of course, nothing is wrong with an 
expansion in powers of the field amplitudes which works when the amplitudes are small but breaks down when they 
are large. On the other hand, it clearly would be desirable to develop a heat kernel expansion in which the Polyakov 
loop enters directly and not only via the generalized zero mode. At the present stage it is not clear whether this is 
possible. 

APPENDIX A: USEFUL SUM-INTEGRALS 

The following sum-integrals show up in the calculation of the lowest order terms of the heat kernel expansion: 

(^rf^-^e-^ =: I {^T) , (Ala) 

^ D ' 2 f ^ 2rfc o =: hti^T) , (Alb) 

^ D ' 2 f£^D e ~ Tk22rk ^ = : MV^)%, (Ale) 

^ D/2 f (Cp e ~ T * r 2fc ° 4 =: ht{ ^ T) ' (Ald) 

(4 ^ )B/ ^(tp e " k2 ±T 2 klhki =: I im (V^T)6 t3 , (Ale) 
{Attt) d / 2 ^ e- Tfe2 AT 2 hk 3 kik m =: h s (s/^T) (SySim + $u6 jm + 5 im 5 ol ) . (Alf) 

Note that sum-integrals with an odd number of k factors vanish. The spatial fc-integrations can be readily evaluated 
using 

f d D k P . . . 

j = 1 ( A2a ) 

d D k Z2 . . 1 



One gets 



■ e k k t kj = ~<%, (A2b) 
d D k -2 1 

— Q-e" fc kikjklkm = - (Sijdlm + SilS jm + SirnSjl) . (A2c) 



Jo(V^T) = hsiV^T) = UtWrT) = {^rf' 2 T £ exp(-r^) , (A3a) 

71 — — OO 

+ oo 

ht(V^T) = hmiV^T) = {AutY^T J2 exp(~TCjl)2rcj 2 n , (A3b) 

Ti — - — OO 

+ oo 

h t (V^T) = (inr^T £ exp(-r<4) -r 2 ^ . (A3c) 
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These finite-temperature sums can be rewritten using a generalized theta function transformation. We adopt formula 
(3.11) from H (with a misprint corrected): 



+00 



2 cxp(- 7 (n + l/2) 2 ) [ 7 (n + l/2) 2 ] P 



n— — 00 



with 



1/2 



(2p-l)! 
2 p 



2 l p! 



l + 2£(-l)" exp(-7r 2 n 2 / 7 )£ 



2 2 \ ^ 
— 7rrr 



p a non-negative integer and (2p — 1)!! := 1 • 3 



l\(p-l)\{2l- 1)!! ' 

• (2p - 1). We apply JHJ| to (JX3J with 7 = (2ttT) 2 t: 



I = l + 2^(-irex P - — 

n=l v 7 

/» = l + jfnrexp^) (1-^). 



/« = l + 2^(-irexpl- i ^) 
n=l v 7 

We also present some combinations of /...-functions which appear in l|35|) : 



1 -4 



4rT 2 3 V4tT 2 / 



4 / n 2 y 



* -26, +4. -J f(-D" exp(-j^) (^) : 

n— 1 v / \ 



Cf 



(A4) 
(A5) 

(A6a) 
(A6b) 
(A6c) 



(A7a) 
(A7b) 



APPENDIX B: PROPERTIES OF THE GENERALIZED ZERO MODE 

In this appendix we collect and prove some important properties for the generalized zero mode defined in l|26|l . 
We first note that for an abelian gauge theory ao is just the zero mode of the Fourier decomposition of the time 
component of the vector potential: 



A (x ,x) = ^2 a n(x)< 



(Bl) 



For a non-abelian gauge theory such an easy relation cannot be obtained. The abelian case, however, justifies our 
name choice "generalized zero mode" . The relations presented in the following hold for arbitrary gauge theories. 
Periodicity: 



a (x + /?, x) — — TlogPexp 



xo+2/3 



dtA (t,x) 



x +P 



-TlogPexp 



xo+P 

dt A (t + /3,x) 

xo =A (t,x) 



a (x ,x). (B2) 



Behavior under gauge transformations: 
We restrict ourselves to gauge transformations which are periodic in xo- These are the only ones appropriate for the 
fermion fields which are integrated out via the heat kernel method. On account of (|20|l we find 

U[A}(x + (3, x; x , a?) — > V(x + /?, x) U[A}{xq + /?, x; xq, x) V~ 1 {x a , x) = V(x , x) U[A](x + f3, x\ x , x) V~ (xo, x) . 

(B3) 

To obtain the transformation property of ao we decompose U as 



U[A](x + 0,x;x o ,&) =: 1 + 



(B4) 
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and use a Taylor expansion of the log-function: 

oo 

logU[A](x +P,£;xo,x) -> logfVW -1 ) = logfy (1 + O)^ 1 ) = log(l + VOV- 1 ) = V (-1)" +1 - (t/OF" 1 )" 

* — ' n 

n—l 

= Vlog(l + 0)U~ 1 = V" (log IT)!/- 1 ■ (B5) 



Thus 

Finally one can also show using the Taylor expansion of the exponential function: 

exp (aao(x)) — ► exp (aV(a;) ao(x) V~ (x)) = V(a;) exp (aao(x)) V~ 1 (x) 
for arbitrary numbers a. 



(B6) 
(B7) 



APPENDIX C: PROOF OF EQUATION (ISHl 

In 0| the derivative of a straight line gauge link Us[A](b,a) with respect to its endpoints a(x), b(x) has been 
considered. One gets 



d*U s [A](b,a) 



db v 
dxu 



A v (b) -(b-a) a J da a U s [A] (b, z) F av (z)U s [A](z,b) 
o 



U s [A}(b,a) 



U 3 [A](b,a) 



A v {a) ~{b-a) a Jds{s-l) U s [A](a, z) F av {z) U s [A]{z, a) 
o 

with the straight line path connecting a and b, 

z = z(a) = a + a (b — a) , 

and the field strength 

F„„ := [D^D V ]. 

This implies for the gauge links of interest 

[D%,U[A}(x o + 0,x;x o ,x)} = O 

and 

D%U[A](x; x') = U[A](x; x')c% . 
We make use of the decomposition of U given in l|B4|) . Using (|C4|I we find 

[D ,O]=0 

and therefore 

[D ,O n }=0. 

This yields 



da u 

OX,, 



oo 1 

[D , a ) = -T [D , log(l + O)] = -T V(-l)" +1 - [D , O n ] = . 



(CI) 

(C2) 
(C3) 
(C4) 
(C5) 
(C6) 
(C7) 

(C8) 



n—l 



Hence we obtain 

oo ^ oo 1 OO oo 1 

[D , exp (x a (x))} = £ - [D , = £ - [D , z«] a£ + £ - afl [A,, ag] = £ - nx^ 1 a£ 



n=0 



a exp(x a (x)) 



(C9) 



Finally we get 



D x U[A](x;x') = D%e Xoa °Mu[A](x;x')e-< ao( - x ^ 

= D x : e x ° ao{x) ] U[A](x;x')e- x '° ao( - x "> + e x ° aa ^ D x U[A](x-x') 
= a (x)U[A](x;x') 

where we have used l|C5l) in the last step. Using (|C10|1 . (|C8j) and (j29a|) leads to 
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